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ABSTRACT 

We present a characterization of the close environment of GRB 980425 based on 5-160 /im spectro- 
imaging obtained with Spitzer. The Gamma-Ray Burst GRB 980425 occurred in a nearby (z = 0.0085) 
SBc-type dwarf galaxy, at a projected distance of 900 pc from an HII region with strong signatures 
of Wolf-Rayet (WR) stars. While this "WR region" produces less than 5% of the £?-band emis- 
sion of the host, we find that it is responsible for 45±10% of the total infrared luminosity, with a 
maximum contribution reaching 75% at 25-30 fim. This atypical property is rarely observed among 
morphologically-relaxed dwarves, suggesting a strong causal link with the GRB event. The lumi- 
nosity of the WR region (is-iooo/jm =4.6 x 10 8 L Q ), the peak of its spectral energy distribution at 
< 100 /im and the presence of highly-ionized emission lines (e.g., [Nelll]) also reveal extremely young 
(<5Myr) star- forming activity, with a typical time-scale of only 47Myr to double the stellar mass 
already built. Finally, the mid-IR over i3-band luminosity ratio in this region is substantially higher 
than in star- forming galaxies with similar Lir, but it is lower than in young dust-enshrouded stellar 
clusters. Considering the modest obscuration measured from the silicate features {rg.^^m ~ 0.015), 
this suggests that the WR region is dominated by one or several star clusters that have either partly 
escaped or cleared out their parent molecular cloud. Combined with the properties characterizing the 
whole population of GRB hosts, our results reinforce the idea that long GRBs mostly happen within 
or in the vicinity of relatively unobscured galactic regions harboring very recent star formation. 
Subject headings: infrared: galaxies — galaxies: individual (GRB980425) 



1. INTRODUCTION 



It is now widely believed that Long Gamma-Ray 
Bursts (LGRBtrM) ori ginate from the cataclismic collapse 
of massive stars~(e.g., Woosley|1993 Galama et al.|1998" 
' - '• "")i)f ' ' 



MacFadyen fe Woos cyj |1999[ |Daigne fc Mochkovitch 
20001 |Heger et al.||2003[ |Stanek et al.||2003| |Meynet fc 
.\ larder 2005) in star-forming regions of distant galaxies 



Fruchtcr 1999 



Bloom et al. 



20021 iHjorthet al.|2002 



Fruchter et al.|2006 )~TEe^ can be observ ed at extremely 



high redshift (z > 6, Greiner et al.|2009[ ) and because of 
their association with young and short-lived stars, they 
are invaluable tools for sign-posting the activity of star 
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formation back to the earliest epochs of cosmic history 
and galaxy evolution. Furthermore, their optical and 
near-infrared afterglows can be used as transient back- 
ground sources for probing the intergalactic space and 
the interstellar medium o f galaxies up to very large cos 



mological d istances (e.g., Vreeswijk et al. 2001 Totani 



et al. 20061. As a result they might eventually provide 



unique constraints on the era of reionization and t he for 
mat ion of the ver y first structures in the Universe ( Lamb 



& Reichart||2000|) 



However, it appears that LGRBs may only trace re- 
gions of star-forming activity in a biased way, in the sense 
that their formation is more likely to occur in low-mass 
and sub-L* galaxies with low-metallicity environments 



( Soko lov et aL"1|20011 |Fynbo et al.||2003| |Courty et al 
| 2004| ,1^0071 [kistler et al.|MQ[ |Kocevski et al.|2M9l |Han 



|et al.|2010[|Niino et al.|201ip . In the latest stages of their 
evolution, massive stars with metal- poor envelopes retain 
a higher angular momentum and they are less subject to 
mass loss. After the final collapse these physical condi- 
tions lead to the formation of fast-rotating black holes 
and the subsequent accret ion of material thus favors the 
formation of GRBs (e.g., |MacFadyen fe Woosley||1999| 



Hirschi et al.||2005[ ). In fact, observational evidence tor 
a possible relationship between LGRBs and chemically- 
young galaxies has already emerged from a number of 
studies. On average the hosts of LGRBs are character- 
ized by blue colors, rather low extinction and modest 
bolometric luminos ities, as well as diffuse and irregular 



morphologies (e.g., Bloom et al. 



2003) |Le Floc'h et al.||2iJ06| |Fruchter et al-pOO^ . TKey 



2002 ] |Le Floc'h et al. 



also exhibit higher Lya emission and larger specific star 
formation rates (SSFR) with respect to other galaxies at 
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simila r redshifts (Fynbo et al. 2003 Christensen et al. 
2004 1 , and their oxygen abundance is lower compared to 
field star -forming sources wi t h similar masses and lumi- 
nosities ( |Modjaz et~aLl|2008[ [Han et al.||2010[ ). Further- 
more, the connection between GRBs and Type Ic super- 
novae has led to the idea that the progenitors of long 
GRBs should more likely or iginate from binary systems 
( Podsiadlowski et al. 2004), which complicates further 
the link that exists between LGRBs and the cosmic his- 
tory of structure formation. Better understanding the 
physical conditions and the properties of the environ- 
ments favoring the formation of these cosmic explosions 
appears, therefore, as a critical step before allowing a 
full exploitation of their hosts as cosmological tracers of 
galaxy evolution. 

Follow-up of Gamma-Ray Bursts has shown that high- 
redshift LGRB hosts are challenging to characterize in 
detail, since t hey tend to be low-m ass and sub-luminous 
systems (e.g., Savaglio et al.[ 2009). Host galaxies at very 
low redshift may thus provide us with some of the best 
cases to explore the intimate connection between the oc- 
curence of GRBs and the properties of their very close en- 
vironment. For example, the now-popular GRB 980425 
was identified at z — 0.0085 from i ts apparent associa - 
tion with the hypernova SN 1998bw (iGalama et al. 1998 ). 
With a luminosity distance of only 36 Mpc it is currently 
the closest GRB known with a confirmed redshift in the 
local Universe. At this distance an angular separation of 
1" on the sky corresponds to 170 pc. It provides there- 
fore an excellent prototype for studying the immediate 
surrounding of LGRBs within their host galaxies. 

The host of GRB 980425 is identified as ESO 184-G82 
in the survey of the sout hern sky that was performed 



by |Holmberg et al. (1977). Observations with the Hub- 
ble Space Telescope (HS'l) showed that it is an isolated 
and barred SBc-type sub-luminous galaxy (Lb = 0.02 L B ) 
harboring a large n umber of active star-forming regions 
( Fynbo et al.||2000[ ) . While the surface brightness of the 
area where GRB 980425 was observed is not particularly 
large, a remarkable result is that the GRB occurred at 
a projected distance of only ~ 900 pc from a quite lu- 
minous HII region characterized by very distinct prop- 
erties compared to the rest of the system. This region 
is the brightest complex of star formation in the host 
and although it only makes a few percent of the total 
luminosity at optical wavelengths, infrared (IR) imag- 
ing obtained with the Spitzer Space Telescope revealed 
that its 5-25 jjm spectral energy distribution (SED) is 
characterized by a steeply-rising contri bution with re 



spect to the whole SED of the galaxy (Le Floc'h et al. 
20061. At 24 /im it is responsible for a fraction as large 



as ~ 75% of the monochromatic luminosity of the en- 
tire system. Among low-luminosity star-forming disks 
and morphologically-evolved dwarf galaxies, such a high 
contribution of extra-nuclear star formation is rarely ob- 
served. It makes the host of GRB 980425 an atypical 
object of the local Universe. 
While the stellar m ass already assembled in t his region 



is low (~ 6 X 10 6 M , Michalowski et al. 2009|), long-slit 



spectroscopy revealed the presence of numerous Wolf- 
Rayet (WR) stars as well as prominent Ha and H/3 emis 



sion lines with very large equivalent widths (Hammer 
ct al. 2006). A mctallicity of Z ~ 0.5 Z was inferred from 



et al.|2006 1 and the estimated star formation rates (SFR) 
vary between 0.06 and 0.3 M ff yr" 1 \he Floc ' h et al.|2006 



Hammer et aL|2006[|Christensen et al.|2008[|Michalowsk'i 
et al.||2009 ). This implies a large SSFR and suggests 
that this environment (denoted as the "WR region" here- 
after) is experiencing a young and vigorous sta rburst. 
Finally, integral fie ld spectroscopy performed by |Chris 



tensen et al. (2008) with the VIMOS instrument revealed 
that the ionizing conditions in this area are markedly 
different from those observed in the other HII regions 
found in the host. With a particularly strong [OIII]/H/3 
line ratio (>5), the ionized gas properties observed in 
the WR region are typical of very young starbursts at 
sub-solar metallicities, and they also appear very similar 
to the properties characterizing the more distant GRB 
host galaxies as a whole. Along with the detection of 
WR stars considered as very strong candidates of col- 
lapsar progenitors, the emission line diagnostics suggest 
therefore an intimate connection between the characteris- 
tics of this source and the physical mechanisms triggering 
the for mation of LGRBs in the Universe. IHammer et al.l 
( 2006 1 even proposed that the WR region may have been 
the true birth place of the GRB 980425 progenitor. The 
latter would have been dynamically ejected from this en- 
vironment, to finally produce a collapsar after traveling 
at velocities ~ 200-300 km s _1 during a few Myrs through 
the interstellar medium (ISM) of the galaxy. 

Following our previous characterization of the 
GRB 980425 host b ased on 4.5 /im 8 /m i and 24 /mi 
broad-band imaging ( Le Floc'h et al.|2006 l we performed 
additional IR follow-up with the Spitzer Space Telescope 
in order to further constrain the spectroscopic mid-IR 
properties as well as the far-IR luminosity of the WR re- 
gion and the host galaxy. Here we report on the 5-35 /im 
spectrum that we obt ained with the In frared Spectro- 



graph for Spitzer (IRS, Houck et al.|2004b[ ) as well as the 
70 /im and 160 /im imaging obtained with the Multi-band 



that 



Imager and Photometer for Spitzer (MIPS, Rieke et al. 
2004|. Along with the host of GRB 0312 03 \z = 0.1055) 
was also observed with the IRS by | Watson et al | 



( |2011| , ESO 184-G82 represents one of the very few iden- 
tifications of LGRB host galaxies located close enough to 
allow detailed characterization in the thermal infrared. 
Our observations are described in Sect. 2, while in Sect. 3 
we present the broad-band photometry and the spectro- 
scopic measurements derived from the Spitzer data. In 
Sect. 4 we derive the total IR SEDs and star formation 
rates of the host and the WR region, and in Sect. 5 we 
constrain in more detail their mid-IR spectroscopic prop- 
erties. We discuss our results in Sect. 6 and we present 
our conclusions in Sect. 7. Throughout this work, physi- 
cal quantities are computed assuming Hq = 70 kms" 1 as 



well a s the Initial Mass Function (IMF) from Salpeter 
(fl955|. 



the effective temperature of the ionized gas (Hammer 



2. OBSERVATIONS AND DATA REDUCTION 

2.1. Far- Infrared imaging 

The host galaxy of GRB 980425 was imaged at 70 /im 
and 160 /im with the MIPS instrument as part of a 
Spitzer Cycle 3 General Observer program (PID: 30 251, 
AOR keys 24966656 & 24966400). The observations were 
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carried out in 2008 using the MIPS "Small Scale Photom- 
etry" mode. With this operating mode a set of 10 frames 
is taken for each observing cycle and a dither is automat- 
ically performed around the targeted source. We respec- 
tively obtained 5 and 12 cycles at 70 /xm and 160 /im, 
which resulted in a total on-source integration time of 
545 s at 70 /im and 250 s at 160 /im. 

For a very modest amount of additional requested time 
we also performed a new observation of the GRB 980425 
host at 24 /im, with the goal of checking if some fraction 
of the luminosity of the bri ght 24 /im point sour ce de- 
tected in our previous data ( |Le Floc'h et al.|[20 06 ) could 
be associated to the emission from a transient object 
(e.g., supernova). A single cycle of 10 frames was ob- 
tained, leading to a total integration time of 165 s on 
source. 

Each individual image was processed with the most 
recent version (v. 3.1) of th e MIPS Data Analysis Tool 
(DAT, [Gordon et at||2005[ ). This update includes ad- 
ditional processing steps with respect to the first ver- 
sion of the DAT initially released. In particular a 
pixel-dependent background subtraction was applied to 
each individual frame before co-adding the images. We 
used calibration factors of 702 (MJysr _1 )/(ADU) and 
41.7 ( MJysr~ 1 )/(ADU) at 70 /im and 160 um r espec- 



tively ( Gordon et al. 



2007 Stansbcrry et al. 2007). The 



pixel size ot the final mosaics was set to halt the pixel 
size of the MIPS detectors and the absolute astrometric 
registration was tied to the Spitzer astrometry of the raw 
data as provided by the Spitzer Science Center. The im- 
ages obtained at the 3 MIPS wavelengths are presented 
in Fig. 1. For reference they are compared with the image 
of the GRB 980425 host galax y obtained at 4.5 /mi from 
our first Spitzer observations ( Le Floc'h et al.|2006 1 . 



2.2. Mid-Infrared spectroscopy 

The mid-IR spectroscopic follow-up of the WR region 
was performed as part of the IRS Guaranteed Time Ob- 
servation (GTO) program on 19 April 2006 (AOR key: 
12246528). The target was first acquired using the red 
peak-up (RPU) camera. During this acquisition an im- 
age of the source at 22 /im was obtained in Double Cor- 
related Sampling (DCS) mode to locate the mid-IR cen- 
troid of the target, and then offset to the appropriate 
slit (see the Spitzer Observers Manual for more details). 
Using the Short-Low (SL, 5.2-14.5 /im) and Long-Low 
(LL, 14.0-38.0 /im) modules of the IRS, we obtained a 
5^38 /im low resolution spectrum of the WR region. The 
total on-source exposure time was 244 s (2 cycles of 60 s) 
per SL order and 126 s (2 cycles of 30 s) per LL order. 
In addition, higher resolution observations were obtained 
with the IRS Short-High (SH, 9.9-19.6 /im) for a total on- 
source time of 975 s (4 cycles of 120 s). For each mode of 
our IRS observations the orientation of the slit aperture 
with respect to the GRB 980425 host galaxy is illustrated 
in Fig. 1 (top-right panel). 

We begun analyzing the observations by processing 
the data with the Spitzer Science Center pipeline (ver- 
sion 15.5) and the extraction of the signal was carried out 
with the Spectral Modeling, Analysis, Reduction Tool 
(SMART Ver. 5.5.1, |Higdon et al |[2004] ). In the case of 
the low-resolution data though, the final reduced spec- 



trum later used in this paper was retrieved from the re- 
cently pub lished Cornell AtlaS of Spitzer IRS Sources 
(CASSIS, |Lebouteiller et al.||2011a| ), which provided us 
with an improved and updated processing of the data. 

In brief, the IRS data reduction starts from the inter- 
mediate pipeline products (droop files), which only lack 
stray light and flat field correction. Individual pointings 
to each nod position of the slit were co-added using me- 
dian averaging. For SL and LL spectra, we considered 
the difference of the two nod positions to remove the con- 
tribution of the background. Then we extracted the spec- 
tra for each nod position using a variable width aperture, 
which scales the extraction aperture with wavelength to 
recover the same fraction of the diffraction limited in- 
strumental Point Spread Function (PSF). The data from 
SH were extracted using the full slit extraction method 
from the median of the combined images. Since no sky 
(off-position) measurements were taken for the high res- 
olution module, the contribution of the sky emission was 
not subtracted from SH spectra. Then the spectra were 
flux calibrated by multiplication with the Relative Spec- 
tral Response Function (RSRF), which was created from 
the IRS standard stars a Lac for SL and LL and £ Dra 
for SH, for which accurate templates were available (Co- 
hen et al. 2003). We built our RSRFs by extracting the 
spectra of the calibration stars in the exact same way as 
the science target, and dividing the stellar templates by 
the extracted stellar spectra. We produced one RSRF for 
each nod position in order to avoid systematic flat field 
errors. The signal difference between the nod positions 
provide the error estimates. Finally, the flux calibrated 
spectra of each order of the low-resolution modules were 
scaled to the first order of LL (LL1, 20-36/iin), which 
was used to define the overall continuum of the source. 

The final low-resolution IRS spectrum of the WR re- 
gion is shown in Figure [2] Given the star-forming nature 
of GRB host galaxies we compared it with t he IRS obser 
vat io ns of the starburst galaxy NGC 7714 ( Brandl et al 



20041 and the super- lu minous HII region IN GC 5461 in 
thespiral galaxy M 101 (|Gordon et al.||2008|). The three 



SEDs exhibit striking similarities but also noticeable dif- 
ferences. First, they are all characterized by a steeply- 
rising continuum of hot dust emission, and the Polycyclic 
Aromatic Hydrocarbons (PAH) features commonly seen 
in star- forming galaxies ( Laurent et al.|2000||Brandl et ah 
2006 Smith et al.|2007 ) are also clearly apparent. These 



PAHs as well as the continuum emission originate from 
dust grains stochastically heated by the radiation field 
and by the UV photons from young populations of mas- 
sive stars. Their luminosity correlates very well with the 
total IR luminosity measured between 8 /im and 1000 /im 
( Brandl et al.|200"6| ) and they are therefore considered as 
a fairly accurate and extinction-free tracer of the activ- 
ity of star formation in galaxies (e.g., Rousscl et al.|200T 



Calzetti et aT|2007| [Diaz-Santos et al |[2008[ ). 
However we note that the 7.7 urn PAH teat 



/im far feature in the 
GRB host is substantially weaker than in NGC 7714 
and NGC 5461, while the main other PAHs detected at 
6.2 /im, 11.3 /im and 12.7/im have roughly the same rel- 
ative strengths in the three spectra. Furthermore, the 
mid-IR spectral slopes are almost identical between the 
WR region and NGC 7714, while the mid-IR spectrum 
of the HII region NGC 5461 seems to rise even more 
rapidly at the longest wavelengths. Finally the IRS spec- 
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Fig. 1. — The host galaxy of GRB 980425 imaged by Spitzer at 4.5 fim (top-left), 24 (im (top-right), 70 /im (bottom-left) and 160 (im 
(bottom-right). The four panels have the same field of view, North is to the top and East is to the left. The position of the previously 
identified Wolf-Rayet (WR) region and the location where the GRB was observed (see Sect. 1) are shown in the 4.5 £im image. The physical 
scale is given in the bottom-left panel while the beam size of the instrument is indicated for each of the MIPS observations. The 24 /im 
image also illustrates the position of the slit apertures for the IRS observations performed with the low (blue dashed lines) and high 
resolution (white dashed line) modules. In each panel the white 'cross' symbol (+) pinpoints the center of mass of the galaxy defined from 
the 4.5 fim image, which shows that the strong emission detected at 70 /im and 160 /jm is shifted toward the WR region. 



trum of the WR region reveals the presence of promi- 
nent fine-structure emission lines such as [SIVJ-10.51 /im, 
[NeIII]-15.55 /xm and [SIII]-18.71 /Ltm, which are also de- 
tected in NGC5461 but either absent or much fainter in 
NGC 7714. In the following sections we will quantita- 
tively analyze the relative strengths of the different fea- 
tures characterizing the WR region of the GRB 980425 
host. In particular we will further explore how it differs 
from what has been so far observed in the mid-IR spectra 
of star-forming galaxies. 

3. MEASURED QUANTITIES 

3.1. MIPS broad-band photometry 

This section describes the broad-band photometry 
measurements obtained from our MIPS imaging. All the 
measured fluxes are summarized in Table[T] 

3.1.1. MIPS-70fim 

As one can see from Figure[T] the host of GRB 980425 
is clearly detected at 70 /im, and in spite of the mod- 
est spatial resolution of the data, low surface brightness 
emission can also be observed up to a distance of ~ 35" 
from the central region of the galaxy. To determine the 
total 70 /im flux density of the host we first estimated the 
sky background in the image by taking the median value 
of the pixels in the outskirt regions up to ~ 100" from 
the center of the source. We then subtracted this sky 
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Rest-frame wavelength (,um) 

Fig. 2. — The IRS low-resolution spectrum of the WR region in 
the host galaxy of GRB 980425 (black solid line), compared with 
the spectral energy distribution of the starburst galaxy NGC 7714 
||Brandl et al.|2004| red solid line) and the SED of the HII region 
NGC 5461 in IVl 101 |Gordon et al.|2008) blue solid curve with thick- 
ness indicating the la uncertainty). Ine light shaded area around 
the GRB host spectrum represents the la error bar estimated for 
each resolution element. Similar to NGC 7714 and NGC 5461 the 
WR region exhibits a steeply rising continuum and bright PAH 
features, but also quite prominent ionic emission lines. 
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TABLE 1 

Infrared photometry and global properties of the WR region and the GRB host 



WR region GRB 980425 host (total) 



Reference 



Fs^rn 
FrOfim 

Fieoiim 

hog(L IR /L e ) 

SFR 



1.8±0.1mjy 
19.5±2.0mjy 
140±25mjy 

8.66±0.04 
O.12M yr- 1 ( a ) 
0.015 



11.9±0.3mJy 
26.2±1.3mJy 
230±30mJy 
615±200 mjy 
9.01±0.07 
0.16 M yr- 1 ^) 



|Le Floc'h et al.| ( |2006| ) 
Sect, d.l.i 
Sect. 3.1.1 
Sect. 3.1.2 
Sect. 4.1 
Sect. 4.2 
Sect. 5 



a Using the calibration from 
b Using the calibration from 



Bell et a.l.|([2005t 



Calzetti et al. (20071 



level to each pixel of the mosaic and we measured the 
photometry inside a circular aperture with a diameter 
of 70" centered on the nucleus of the host. The photo- 
metric uncertainty was estimated from the pixel-to-pixel 
variations measured over the sky. We derived a total flux 
of S 70fim = 230±30 mjy. 

Furthermore, the peak of the 70 /mr emission is clearly 
shifted from the center of mass of the host defined from 
the 4.5 /im stellar emission, and there is strong evidence 
that the bulk of the 70 /im luminosity originates from the 
WR region that also dominates the mid-IR luminosity 
of the galaxy. Figure[3] shows a larger field of view of 
the MIPS-24 /im image of the GRB 980425 host, overlaid 
with the contours of the emission detected at 70 /im. We 
see that the peak of the 70 /im emission coincides with the 
very bright 24 /im point source lying in the South- West 
spiral arm at 12" from the center of the host galaxy. 

Although this distance is slightly smaller than the 
FWHM of the MIPS-70/im PSF (18"), we believe that 
the significance of this shift is robust given the high sig- 
nal to noise of the 70 /im emission detected toward the 
GRB host. As described in the next paragraph the peak 
of the 70 /im emission originates from a point source de- 
tected with a signal to noise S/N ~ 5.6. The positional 
la uncertainty associated with this point source can be 
estimated as 0.6 x FWHM x (S /N)- 1 = 6.75" (see Ap- 
pendix B of Ivison et al.||2007| ), which is twice smaller 
than the measured offset. Also, such a shift with respect 
to the nucleus of the galaxy could not be artificially pro- 
duced by a misalignment between the 24 /im and 70 /im 
mosaics given the astrometric calibration of the Spitzer 
data (la ~ 1" rms). In fact there are two other galaxies 
detected to the South of the GRB host (see Figure[3]) and 
their 70 /im emission appears to be very well centered on 
the peak of their 24 /im counterpart. 

We compared the two-dimensional profile of the 70 /im 
detection of the GRB 980425 host with the MIPS-70 fim 
PSF. In spite of the lower spatial resolution of the MIPS 
data at this wavelength we found that there is a spatially- 
resolved low surface bightness component extending to 
the East of the peak of the 70 /im emission. This com- 
ponent is observed over the whole extent of the galaxy 
and it is clearly apparent from the asymmetrical shape of 
the contours displayed in Figure[3] (note however that the 
small extension to the South comes from a data reduc- 
tion artifact produced along one of the two dimensions of 
the 70 /im pixel array). We thus decomposed the 70 /im 
image of the GRB 980425 into a point source centered 
on the peak of the 70 /im emission and the contribution 
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Fig. 3. — The surrounding of the GRB 980425 host galaxy imaged 
at 24 /im, superimposed with the contours of the 70 fim emission. 
The GRB hos t is visible to the N orth while Galaxy A and Galaxy B 
identified by | Foley et al.| ( |2006[ > are also clearly detected at both 
wavelengths to the boutn Jiast. The contours scale linearly from 
1.7 to 15MJysr — 1 . As in Figure 111 the white 'cross' symbol (+) 
indicates the center of mass of the GRB host galaxy. In spite of the 
large FWHM of the far-IR MIPS PSF the peak of the 70 fim emis- 
sion detected toward the host of GRB 980425 clearly coincides with 
the bright 24 (im point source that corresponds to the WR region. 

from the spatially-extended regions u sing the PSF fi tting 
routine of the DAOPHOT package ( |Stetson||T987| . We 
found that the point source contributes to 140±25rnJy, 
which shows that a fraction as large as ~ 60±10% of the 
70 /im emission of the GRB 980425 host originates from 
the WR region. 



3.1.2. MIPS-160fim 

At 160 /tm the GRB 980425 host is also clearly detected 
in our data (Figure [I}. However the FWHM of the MIPS 
PSF at these wavelengths is comparable to the angu- 
lar extent of the galaxy on the sky (~40"), and given 
the modest signal to noise of the observations it was not 
possible to confirm the presence of a spatially-extended 
component similar to what we found in the 70 /tm im- 
age. The total flux density was therefore estimated with 



G 
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aperture photometry in a radius of 40", and a correction 
taking account of the PSF emission outside of this aper- 
ture was applied f o llowin g the prescriptions described by 



Stansberry et al 



(|2007j). We obtained a total flux of 
i60pm = 615±20U mJy. Given the relatively small size 
of the 160 /im pixel array as well as the contamination 
from the galaxy detected to the South of t he GRB host 

-if- 

of 



tion as well as on the strength and the hardness of their 
radiation field. To constrain the luminosity of these fea- 
tures we analyzed our low-re solution spectrum w ith the 



(referred as Galaxy A by |Foley et al. 
: [3]), only a very small 



2006 see also Fig- 



pixels were avail- 



ure only a very small number 

able to estimate the sky background in the image. This 
explains the large photometric uncertainty affecting our 
flux estimate. 

Given the low angular resolution at these long wave- 
lengths it was not possible to perform a spatial decom- 
position similar to the one described above for the 70 /im 
data. In particular, we could not determine the frac- 
tion of energy powered by the point source that domi- 
nates the monochromatic luminosity at shorter IR wave- 
lengths. However, FigureO] reveals that the peak of the 
160 /im emission is again shifted to the South- West with 
respect to the center of the galaxy. Within the astro- 
metric uncertainties it coincides with the peak of emis- 
sion detected at 24 /im and 70 /im, which suggests that 
the luminous HII region detected in the mid-IR is also 
responsible for a very large fraction of the total far-IR 
emission of the galaxjj^] 

3.1.3. MIPS-24fim 

The new 24 /im image of the GRB 980425 host galaxy 
obtained in 2008 as part of this program looks very sim- 
ilar t o the observations pe rformed by MIPS / Spitzer in 
2004 ( |Le Floc'h et al.||2006|) . Most of the emission origi- 
nates from a single point-like source that appears to co- 
incide with the WR region, while a couple of other HII 
regions as well as a more diffuse component are also de- 
tected toward the nucleus and in the spiral arms of the 
host. We measured the total flux of the galaxy and the 
luminosity of the WR region using a PSF fitting decom- 
position and the same approach followed for the 70 /im 
data. We found that the unresolved component con- 
tributes to a level of 19.5±2.0mJy, while the total emis- 
sion of the galaxy reaches = 26.2±1.3 mJy. Within 
the error bars these r esults are fully consiste nt with the 
24 /im photometry of |Le Floc'h et al.| ( |2006l ). The lack 
of temporal variations shows that if some dust heating 
transient emission produced after the GRB explosion has 
contributed to the mid-IR emission of the galaxy, this 
emission must either be very faint (i.e., <10% of the to- 
tal luminosity) or decrease very slowly with time. 



3.2. Mid-Infrared spectroscopic features 

Numerous emission lines are clearly visible in the 
IRS spectroscopic observations of the WR region (Fig- 
ure^. They are commonly observed in the spectral en- 
ergy distribution of star-forming galaxies and they pro- 
vide valuable information on their activity of star forma- 

12 The host of GRB 980425 is the only object covered by the field 
of view of the 160 /im observations. Therefore it was not possible 
to independently confirm the astrometry with other field galaxies 
as we did at 70 (im and we had to rely on the Spitzer astrometric 
registration which is accurate to 1". 



PAHFIT tool developed by Smith et al. (20071. This 



code was specifically designed for decomposing Spitzer 
IRS spectra into the contribution of the different com- 
ponents that characterize galaxy mid-IR SEDs. These 
are the hot dust continuum produced by the Very Small 
Grains (VSG), the broad-band PAH features, the narrow 
ionic emission lines and the effect of silicate absorption 
at 9.7 /im and 18 /im. The result of our decomposition is 
displayed in Figure [4] Given the large number of indi- 
vidual features identified with PAHFIT, only the spectral 
profile corresponding to the main PAHs is represented, 
along with the shape of the underlying continuum (thick 
solid line). Similarly, only the position of the most promi- 
nent forbidden lines is labeled on the plot, and the total 
SED reconstructed from the model is shown as a blue 
solid line. 



100.0 



10.0 7 




Rest-frame wavelength (/im) 

Fig. 4. — Decomposition of the IRS spectrum observed in the 
WR region o f the G RB host, obtained with the PAHFIT tool of 
|Smith et H] ( |2007| l. The measured SED is illustrated with the 
red data points and the associated error bars, while the total SED 
reconstructed from the model is shown with the blue solid line. The 
underlying continuum and the modeling of the main PAH features 
at 5.3, 6.2, 7.7, 8.6, 11.3, 12.7, 14.2 and 17 /mi are represented 
by the black solid curves. The position of the different ionic lines 
identified by the fit are also indicated. 



The narrow ionic lines observed between 10 /im and 
19 /im are also clearly visible in the high-resolution spec- 
trum that we obtained with IRS. A more accurate esti- 
mate of their luminosity was therefore obtained by fit- 
ting a gaussian function to their emission profile derived 
from the SH data. This is illustrated in Figure 5j which 
shows the SH IRS observations obtained for the [SIV], 
[Nell], [Nelll] and [SIII] ionic lines. Since no sky sub- 
traction could be performed for these data though, the 
underlying continuum and the equivalent widths (EW) 
of these features were estimated from the continuum of 
the low-resolution spectrum as modeled with PAHFIT. 
Our measurements are given in Table [I] along with the 
fluxes and the equivalent widths of the main PAHs and 
ionic lines measured in the low-resolution data. For the 
strongest and isolated PAHs (i.e., 6.2 /im, 7.7 /im, 8.6 /im 
and 11.3 /im), we indicate the results obtained with the 
global PAHFIT decomposition but we also provide the 
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measures that we derived with a local fit of the contin- 
uum underlying each individual feature using a spline 
function. The latter approach has been commonly used 
in the literature to characterize the mid-IR spectra of 
star-forming galaxies. It usually leads to lower values 
than obtained with PAHFIT, since PAHFIT accounts 
for the full extent of the PAH wings. 



TABLE 2 
Mid-IR Emission Features 
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Fig. 5. — A close-up view of the profiles of the [SIV], [Nell], 
[Nelll] and [SIII] ionic lines observed with the IRS Short-High mod- 
ule. In each panel the observed data corrected for the redshift of the 
source are shown with the black tickcd-solid curve while the model 
of the emission profile is performed with a gaussian fit illustrated 
in red; the expected position of the forbidden line is indicated by 
the vertical dashed line. 



4. TOTAL IR SEDS AND STAR FORMATION RATES 

4.1. Infrared SED fitting 

We constrained the spectral energy distributions of 
the GRB host galaxy and the WR region over the full 
infrared wavelength range by fitting the Spitzer broad- 
band photometry with the empirical lib r aries o f galaxy 
templates pu blish ed by|Chary fc Elbaz (2001), Dale fc 
Helou (2002) and Lagache et al. 12004), as well as with 
the phy sical SEDs derived from radiativ e transfer mod- 
eling by Siebenmorgen & Kriigel (20071. Between these 
different libraries the SEDs mostly vary in the relative 
strength of the PAH features with respect to the hot 
dust continuum, as well as in the temperature and the 
emissivity of the cold dust component sh aping the peak 
of the SEP in th e far-IR. In the library of Siebenmorgen 
& Kriigel (2007) the SEDs also depend on the size of 
the star-forming region responsible for the IR emission, 
and they are given for radii of 0.35, 1, 3, 9 and 15kpc. 
We only considered sizes of 0.35 kpc and 3kpc for the 
WR region and the whole galaxy, respectively. 
To obtain the best possible constraints we combined 



the MIPS fluxes presented in Sect .JO with the IRAC 
8 ixm photometry already published by 
(2006). The fitting was performed separately for each 



Feature 


A 


Flux 


rliVV 




(Cm) 


X 1U (W cm J 


(Mm) 


PAH a 


6.2 


3.14±0.33 (1.7±0.07) 


0.98 (0.38) 


PAH a 


7.7 


9.89±1.39 (3.7±0.5) 


2.12 (0.43) 


PAH a 


8.6 


1.98±0.23 (0.56±0.2) 


0.34 (0.06) 


PAH a 


11.3 


4.03±0.20 (1.7±0.3) 


0.49 (0.15) 


PAH 6 


12.7 


1.16±0.23 


0.12 


PAH 6 


17 


7.58±0.75 


0.58 


[Aril] 6 


7.01 


0.18±0.16 


0.05 


[Arlll] 6 


9.01 


0.19±0.08 


0.04 


[SIV] C 


10.51 


2.01±0.18 


0.40 


H2(S2) C 
[NeII] c 


12.33 


1.35±0.07 


0.22 


12.82 


1.48±0.11 


0.23 


[NeIII] c 


15.55 


2.68±0.14 


0.34 


[SIII] C 


18.71 


1.57±0.12 


0.16 


[SIII] 6 


33.53 


1.90±0.09 


0.17 


[Sill] 6 


34.86 


1.44±0.09 


0.13 



a Main values result from PAHFIT measurements in the low 
resolution data. Values in parenthesis give another estimate 
independently obtained using a local spline-type fitting of the 
continuum. 

b Measured with the IRS low resolution module using PAH- 
FIT. 

c Measured with the IRS high resolution module, except for 
the continuum which was derived from the low resolution 
data. 



library, using th e code LePhare ( Arnouts et al.|[l999| pi 



bert et al.|20 06). Although most ot the IR SED templates 
from the aforementioned librari es vary as a function of 
the total IR luminosity (but see Dale & Helou 2002 for 
a dependence on dust temperature], their normalization 
was kept as a free parameter and the best templates were 
derived from a basic % 2 minimization of the fit. In the 
case of the WR region, we did not include the photom- 
etry at 160 /im since at this wavelength we were unable 
to separa te its c ontribution from the emission of the host 
(see Sect. 3.1.2). We checked however that the best fits 
obtained ior the WR region did not exceed the total flux 
measured for the GRB host galaxy at 160 /im. 

The results are illustrated in Figure|6 which shows 
the measurements from our broad-band photometry and 
IRS spectroscopy together with the global range of 
possible SED fits defined from the best-fit templates 
that were obtained for each of the 4 libraries. As ex- 
pected from the different broad-band fluxes measured 
with IRAC and MIPS we note that the mid- to far- 
IR spectral slope of the WR region is much steeper 
than observed for the whole GRB host. By integrating 
the best fit SEDs between 8 /im and 1000 /im, we de- 
rived total IR luminosities of log(LiR/L0) = 8.66±0.04 
and log(L IR /L ) = 9.01±0.07 for the WR region and the 
whole GRB host galaxy, respectively. In these estimates 
the uncertainties were obtained by combining (i) the dis- 
persion measured between the 4 best fits obtained from 
the different libraries and (ii) each of the dispersions 
measured within the individual libraries themselves. Our 
results thus suggest that 45±10% of the total IR lumi- 
nosity of the GRB 980425 host galaxy originates from the 
WR region close to which the GRB was observed. 

The lack of photometry at submillimeter wavelengths 
prevents us from fully constraining the Rayleigh- Jeans 
emission tail, which traces the properties of the cold 
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Fig. 6. — The Spitzer broad-band photometry of the GRB 980425 host galaxy (red filled diamonds) along with the IRAC /MIPS photometry 
and the IRS spectrum of the WR region (red filled squares, thick black line). Photometric uncertainties are shown with vertical error bars 
when they are larger than the symbol size. The shaded regions represent the range of possible SEDs fitting the 8-160 fim broad-band 
photometry of the host galaxy and the WR region (blue and green shaded areas, r espect i vely), defined f r om th e 4 best-fit templates 
independently obtained using the libraries of [ Chary fc Elbazl 12001) , [Dale fc Heloul {[ 2002) , |Lagache et al.| <|2004b and JSiebenmorgen fc| 
Krugel (2007) (see Sect. 4.1 for more details). For comparison "the IR SEP ot iNCC 7714 taken trom |Marsh"aIret al.|{2007| l and scaled down 
by a factor of 100 is also illustrated (thin red line). 

dust component in the GRB host galaxy and the WR re- 
gion. Therefore the position of the IR SED peak and the 
characteristic temperature of the dust grain emission can 
not unambiguously be determined with the present data. 
However, since we used different libraries of templates 
that significantly vary from one another, the results of 
our SED fitting suggest that we can reasonably predict 
the entire IR spectral energy distribution of both the 
GRB host and the WR region within an uncertainty of 
a factor of ~ 2-3. This can be explained by the fact that 
the IRAC/MIPS constraints imposed on the Wien side of 
the SED substantially reduced the number of possible fits 
from those libraries. In particular, we note that the SED 
of the WR region peaks at much shorter wavelengths 
than the SED of the wh ole galaxy (as also predicted by 
|Micha lowski et al. 2009 from stellar population synthesis 
modeling), showing that the characteristic temperature 
in the WR region is about twice larger than the aver- 
age dust temperature in the interstellar medium of the 
galaxy. This is also a direct consequence of the larger 
fraction contributed by the WR region at 24 /um than at 
70 /im. 



4.2. SFR determination 

Based on the total IR luminosities computed in the 
previous section we obtained a revised estimate of the 
star formation rates in the WR regi on and the GRB 
host using the calibration provided by Bell et ak| ([2005 ) : 
SFR (M yr" 1 ) = 1.8x lO" 10 (L IR + 3.3L 2800 ) / L . In 
this relation, L2800 represents the monochromatic lumi- 
nosity vL v computed at 2800A, which we derived from a 
linear interpo lation between the NUV a nd U-band fluxes 
published by Michalowski et al. (2009). In the case of 
the WR region we obtained SFR~ 0.12 M yr" 1 . As- 
suming the stellar mass M » — 5.7 x 10 6 M inferred by 



Michalowski et al. (20091, our estimate thus leads to a 
specific star formation rate SSFR — 21 Gyr -1 , which cor- 
responds to a typical time scale of ~ 47 Myr for this star- 
forming complex to double its mass. 

Our new estimate of the SFR in the WR region is 
smaller than the one initially published by |Le Floc'h] 
et al. (2006) from the first Spitzer observations of the 
GRB 980425 host (SFR- 0.35 M yr" 1 ), but this first 
determination had relied on a large extrapolation of the 
24 /im flux to infer the total IR luminosity characterizing 
this environment. Our revised SFR is in good agree- 
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ment with that obtained by Michalowski et al. (20091 
from broad-band SED fitting (SFR = 0.10 M yr ijTand 
it is also quite consistent with the values that can be 
inferred from the direct conversions recently established 
between the SFR and the mid-IR emission of galaxies 



(lAlonso-Herrero et al. 2006 Calzetti et al 



Santos et al. 



For example, Alonso-1 



2007 Diaz 



errero et al 



] 2U06| obtained a tight relationship between the 24 /im 
and the extinction-corrected Pa a luminosities of HII re- 
gions and star-forming galaxies in the local Universe 
(SFR /Mm yr' " 



8.45 x lO" 38 [L^/ergs- 1 ] ' 871 ), 



while Calzetti et al. (20071 established a comparable 
correlation between the SIR derived from the Ha re- 
combination line and the 24 /Ltm emission of star-forming 
regions in nearby extragalactic source^ We applied 
these two relations to the 24 /jm flux of the WR re - 
gion measured from our PSF fitting (see Sect. 3.1.3), 
which respectively led to star formation rate estimates 
of SFR ~ 0.14 M© yr- 1 and SFR- 0.13 M Q yr" 1 . Con- 
sidering the possible variations of the two calibrations 
with the met allicity and the age o f the stellar popula- 



tions (~20%, |Calzetti et al. 2007), we find therefore a 



very good agreement with the revised estimate reported 
above from the UV and IR luminosities. 
In the case of the en t ire ho st gal axy, the relations from 



Alonso-Herrero et al. 



(2006) and Calzetti et al. 
l 



(2007) 
r 



lar populations ( Savaglio et al.|2009 ) and it is likely that 
a large fraction of the total U V and IR luminosities of 
the host is not directly produced by on-going or recent 
star formation. Conseque ntly, the es t imate of the SFR 
based on the relation from Bell et al. ( 2005 ) is probably 
overestimated. 



5. MID-INFRARED SPECTRAL PROPERTIES 



In Sect.[372]we used the PAHFIT tool to analyze the 
different components contributing to the IRS spectrum 
of the WR region (see Figure|4]). Focussing first on the 
shape of the underlying VSG continuum, the decomposi- 
tion reveals that there is almost no absorption from the 
silicate features at 9.7 fim {r^.j^m ~ 0.015). A ssuming 
A v I A oi.7 llm ~ 20 as observed in the Milky Way (|Mathis| 



19901, the optical depth toward the silicates correspond 
to a small extinction Ay — 0.3 mag. This is actually 
consistent with the modest obscuration that was derived 
toward the very center of the WR region based on the 
Hydrogen Balmer decrement m easured with IFU spec - 
troscopy (E(B-V) ~ 0.03 mag, IChristensen etal | 2008). 
It suggests that most of the activity in this area originates 
from a compact and optically-thin environment, and that 
the contribution from dust-enshrouded star formation 
must be very modest. Note that our estimate appears 



13 Note that the relation inferred by |Calzetti et aL 1 (|2007) needs 
to be corrected by a multiplicative factor ot 1.5t) to match the 
Salpeter IMF assum ed in our analysis. See the Appendix A2 of 
|Calzetti et~aT| (2007) for more details. 



smaller than the value obtained by Hammer et al. ( 2006 1 
using optical long-slit spectroscopy (Ay ~ 175 mag). 
Such differences could origina te from aperture effects 



since the extinction inferred by Christensen et al. ( 2008 1 
reaches E(B — V) — 0.17 mag when it is measured over 
the entire WR region. It could also be due to the gen- 
eral difficulty in properly defining the underlying mid- 
IR spectral continuum around the 9.7/^m silicate feature 
when prominent 7.7 fim and 11.3 /im PAHs are detected. 
This degeneracy in the mid-IR SED fitting results in ad- 
ditional uncertainties affecting the determination of the 
obscuration from the silicates, and it could have led to a 
small underestimate of the extinction in the case of the 
WR region. 

Furthermore, the detection of prominent forbidden 
lines with high ionizing potential such as [Nelll], [SIV] 
or [Arlll] indicates the presence of a particularly hard 
radiation field. We quantified the ionizing conditions 
in the WR region using the lines fluxes measured from 
both the low and high-resolution IRS spectra, which led 
to [NeIII]/[NeII], [Arlll] / [Aril] and [SIV]/[SIII] ratios 
of 1.8, 1.05 and 1.3 respectively. Taking into a ccount 
the sub-solar metallicity (Z<O.5Z ) inferred by |Ham-"| 



mer et al. (2006) and Christensen et al. (2008), these 



give a total SFR of 0.18 Meyr" 1 and 0.16 M fl yr^re- pact 
spectively, while the relation obtained by |Bell et al.| 
( ]2005[ ) leads to a total SFR - 0.66 M yr" 1 . At opti- 
cal wavelengths though, the global SED of the galaxy is 
dominated by t he contribution of pas sively-evolving stel- 



ratios are comparable to what has been observed in 
other individually-resolved HII regions and Blue C om- 
Dwarves with similar chemical abundanc es (Wu 



et al.|| 2006| |Gordon et aL|2008|IHao et al.l|2009[ ). How^ 
lese values are substantially higher than what has 
been typically measured in star-forming and starburs t 



galaxies (e.g., Brandl et al. 2006 Beirao et al. 2008) 



where the emission is likely averaged over different stel- 
lar populations spanning a wide range of properties and 
ages. Assuming that the mid-IR emission of the WR re- 
gion originates from the c entral star cluster that produces 
most of the optical light ( |Hammer et al.|[2006 ), the spa- 
tial extent of the B-band emission (i.e., ^100pc) and 
the hardness of the radiation field derived with IRS sug- 
gest that most of the mid-IR luminosity is pro duced b 
very young massive star formation. Following |Rigby 



Rieke ( |2004[ ) who studied the evolution of various mid-IR 
line ratios with time using the CLOUDY photoionization 
model, we find that the [NeIII]/[NeII], [ArIII]/[ArII] and 
[SIV]/[SIII] ionizing indices measured in the GRB 980425 
host are compatible with a starburst episode younger 
than ~ 5 Myr. This is consiste nt with the estimate ob- 
tained by |Hammer et al.| (|2006| ) based on the number of 
Wolf-Rayet stars observed in this region (i.e., 1-6 Myr). 

The mid-IR emission line diagnostics that can be used 
to characterize the ionizing conditions of star-forming 
environments are obviously not independent from each 
other. A number of correlations have been established 
between [NeIII]/[NeII], [ArIII]/[ArII] and [SIV] / [SIII] , 
as well as with [SIV] /[Nell] and with optical line ratios 
like [OIII1A5007/[OII1A3727 (e.g.JGiveonet al |2002l[Wu 



ei"l^[2006l [Dale et ai.l(2006l Groves et al.||^00§l IGordon 
et al.||2008l |Hao et al.[|2009| Bernard- Salas et al.| |2009 ) . 



The mid-lR spectral properties ot the WR region are 
fully consistent with these different relationships, as well 
as with the rough correlation that exists between the 
metallicity of galaxies and the hardness of the ir radia- 



tion field (Wu et al. 2006 Gordon et al. 2008). More 



over, the electron density m the WR region can be con 



10 
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strained from the two Sulfur III emission lines detected 
with IRS. The flux ratio [SIIL] A1 8 71/[SIIIl A33.5 3 ~ 0.8 



3 (|Giveon et al 

that 



corresponds to N e ~ 200 cm 
This is slightly larger than the estimate 
derived from the optical spectrum of the 
gion using the [S III] A6716A/ [SIII] A673 1 A 



2002) 



was 
WR re- 
flux ratio 



(N e = 158 cm , Hammer et al. 2006) and it is also 



higher than the densities measured in t he HII regions 



of M 101 (SIII] A18.71/[SIII]A33.53~ 0.2, Gordon et al 



20081) . This could be related to the effect ot the stellar 
winds produced by the Wolf-Rayet populations, which 
may have efficiently compressed the gas and the ISM of 
the WR region in the host of GRB 980425. 

Finally, the equivalent widths of the PAH features 
characterizing the WR region fall within expectations 
given the metallicity and the hardness of the radia- 
tion field previously measured. In star-forming galax- 
ies the relative strength of PAHs over the VSG contin- 
uum is known to remain constant over a wide range 
of environments, but it rapidly decreases with harder 
ionizing radiations. This is particularly apparent when 
the hardness of the radiation field reaches a threshold 
of [NeIII1/[N eII1 



Lcbouteiller 



1-3 ( [Gordon et al."||2008[ 
et al.pOllb ), or similarly when the chemical abundance 
gets lower than 12 + log[Q/H| ~ 8-8.3 (e.g ., |Houck et al 



2004a[ |T^ elbrachret~al.||2005| |Wu et aL||2(Mj| |Engel- 
bracht et al.||2008[ ). The properties of the WR region in 
the GRB 980425 host galaxy correspond to this regime 
of transition between extreme metal-poor regions and 
environments with solar metallicity, and the PAHs are 
clearly detected albeit with smaller EWs than observed 
in sources with softer radiation fields. Looking at Fig- 
ure^ we also note that the 7.7/zm/11.3/zin PAH inter- 
band ratio is substantially smaller than observed in the 
mid-IR spectrum of the prototypical starburst galaxy 
NGC 7714 but it is once again in agreement with the 
trend that has been seen between this quantity and the 
ionizing index meas ured with the [Nelll] / [Nell] line ratio 
( |Smith et aL|2007[ ). Given that the 11.3 firm feature can 
be produced by larger PAHs than the one s responsible 
for the 7.7 /im emission ( Draine fc Li||2007 ), the WR re- 
gion could thus be experiencing a selective destruction of 
the small PAHs th at emit at short mid-IR wavelengths 
( [Smith et al]|2007| ). 



6. DISCUSSION 

6.1. Comparison with other known extragalactic sources 

The total IR luminosity derived in Sect. 4 
(log [Lir/L q ] ~ 8.7) places the WR region at the 
very bright end of the luminosity f unction of HII regions 



observed in the local Universe (Bradley et al. 2006 



Lee et al 
optical lig 

the h ost (<100pc, |Fynbo et al.||2000| |Hammer et al 



2011]) . Considering the physical size of the 
it distribu tion measured in the HST image of 



2006 1 , this environment is also one of the single isolated 



Hll regions with the highest star formation density 
identified so far. For instance, the brightest star-forming 
compiex in the Milky Way, W49A, releases up to 
Ltr ~ 2.7 x 10 7 L m within a size less than ~50pc in 



diameter (Ward-Thompson & Robson 1990), while the 
total luminosity of the 30Doradus nebula in th e Large 
Magelanic Cloud only reaches Lir ~4x 10 7 L (Werner 



et~aT1[l978] ). Yet, some of the intrinsic characteristics of 
this environment have already been observed in other 
extragalactic sources, which suggests that such complex 
of star formation does not represent a unique object in 
the local Universe. For instance, isolated HII regions 
with similar mid-IR spectral slope and ionization indices 
have been found in the giant spiral galaxy M 101 



( |Gordon et al.||2008 ). Among these regions, NGC 5461 
does contai n a population of Wolf-Rayet stars ( |Schaerer 
et al.|[l999 ) and it has a 24 /im flux of ~ 920 mJy, which 
given the distance of M 101 (6.7 Mpc) corresponds to a 
mid-IR luminosity 50% larger than what we measured 
in the host of GRB 980425. The mid-IR properties 
of the WR region also resemble what is seen in Blue 
Compact Dwarf (BCDs) galaxies with comparable 
oxygen abundance, which often exhibit a steeply rising 
hot dust continuum with prominent ionic lines and very 
little absorption by the silicate features ( Wu et aL| 2006). 
Similarly, some of the giant Super Star Clusters in the 
overlap region of the Antennae exhibit comparable or 
steeper mid-IR continuum emission as well as higher 
star formation rates, although these clusters may be 
spatially more extended and their radiation field is 
typically a bit soft er than observed in the WR region 
(Brand! et aL][2009 |. 

What is nonetheless very unusual in the host of 
GRB 980425 is the large fraction (45±10%) of the to- 
tal Infrared luminosity of the galaxy contributed by 
the WR region alone. This large contribution to the 
global star-forming activity of the host is not apparent 
at optical wavelengths, because the ratio of the mid- 
IR and optical luminosities in the WR region is much 
higher than measured over the whole spatial extent of 
the galaxy. This is illustrated in Figure[7J which shows 
the 24 /im over i?-band monochromatic luminosity ra- 
tio L MIR /L B = {vL u ) 2iiim l{vL u ) B of the host and the 
WR region, plotted as a function of their total Infrared 
luminosity. In both cases the J 3 -band luminosities were 

For comparison 



(2009) 



Dale et al | ( |2006| ) and |Smith et al. (12007) ii 
Infrared Nearby G alaxies Survey ( SINUS ) , 



taken from Michalowski et al. 

we over plotted the ratios der i ved fr om the Sloan Digital 
Sky Survey by Hwang et al.| fl2010| ) , thos e measured by 

in the Spitzer 
and the ra- 
tios characterizing the thr ee BCD prototype s SBS 0335- 
052, IIZw40 and He2-10 ( |Hunt et al.|j2005[ ). While the 
Lmir/Lb integrated over the host of GRB 980425 is typ- 
ical of quiescent star-forming galaxies with similar lumi- 
nosities, we see that the ratio in the WR region is more 
than an order of magnitude higher, and comparable to 
that observed in Blue Compact Dwarves. 

This characteristic of the WR region is actually con- 
sistent with what can be observed in recent starbursts, 
despite the apparently modest extinction derived from 
the 9.7/xm silicate features and the Hydrogen Balmcr 
decrement (see Sect.Js]). In fact, very young and dust- 
embedded stellar clusters with much higher extinction 
(e.g., Ay > 5 mag) can harbor Lmir/Lb ratios even 
larger than found in the WR region. As an example 
we included in Figure[7| the case of the "hype rnebula" 
starburst enshrouded in the core of NGC 5253 (Alonso- 



Herrero et al.||2004 ) as we ll as the massive star clusters 
identified in NGC 1365 bylGalliano et all p008l. These 



complexes of star formation share strong similarities with 
the WR region in the GRB 980425 host in terms of SFR, 
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age and luminosities. Yet, the optical depth toward these 
sources is much more significant, explaining their ex- 
treme optical-to-infrared colors. Assuming that the ob- 
scuration and the IR/optical luminosity ratio decrease 
with the starburst age, what we are seeing in the WR re- 
gion could thus be one or several star clusters that have 
already evolved and partly escaped the original molec- 
ular cloud where they were born. Alternatively, and 
rather than an "aging" effect, the difference of visual ex- 
tinction between the WR region of the GRB 980425 host 
and the enshrouded clusters of NGC 5253 and NGC 1365 
could also be a direct consequence of the presence of the 
WR stars themselves, which may have accelerated the 
process of clearing out the local ISM of the star-forming 
region via their emission of long-lasting supersonic winds. 
Depending on geometry effects such a scenario could lead 
to small silicate absorption as observed, while a still im- 
portant emission of very hot dust at mid-IR wavelengths 
could arise from amorphous carbon grains heated to very 
high temperatures by the massive stars dominating the 
spectrum of the WR region. 



fl I M I 1 | [ I 1 I [ I 1 1 I | I 1 I I I 1 I II | 1 I I I 1 I II 1 ] II I 1 II I 1 I | I 
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Fig. 7. — The mid-IR over B-band monochromatic luminosity 
ratio (Lmir/Lb) plotted against the total IR luminosity of the 
WR region (red filled diamond) and the entire GRB 980425 host 
galaxy (blue filled square). For comparison we also indic ate the 
local star-forming galaxy samples of SINGS (open square s. |Smith| 
|et al.|2007| and SDSS ('+' symbols, I Hwang et al.|2010j , the pro- 
totypical Blue Com pact Dwarves SBs 0330-052, 11 kw 4U and He2- 
fO (pink triangles, IHu nt et al. 2005), the embedded starburst in 
NGC 5253 ('*' symbol. lAlonso-fierrero et al . 2004) an d the young 
star clusters of NGC 13bo (lower limits, |Ualliano et al.|2008| l. Mid- 
IR luminosities were computed based on Mli- , S-^4 /im fluxes except 
for the SDSS sample where the IRAS-25 (im bandpass was used as 
a proxy. The Lmir/Lb ratio in the WR region is much higher 
than typically found in star-forming galaxies with comparable lu- 
minosities, but it is lower than observed in young star clusters still 
embedded in clouds of dust and molecular gas. 



6.2. Long Gamma-Ray Bursts and super-luminous 
HII regions 

As demonstrated in Sect. 4, the most striking result of 
our analysis resides in the bolometric luminosity of the 



GRB 980425 host, half of which originates from a single 
and compact star-forming complex located at more than 
2 kpc from the galaxy center. This large contribution 
of extra-nuclear star formation has also been observed 
in a few other spiral galaxies of the local Universe (e.g., 



NGC 5257, NGC 6670, Arp256, Haan et al. 2011), al 



though these cases are rare and they have much higher 
mass and SFR than the host of GRB 980425. Most of- 
ten, the star-forming activity of morphologically-evolved 
systems in the luminosity range of LGRB hosts appears 
to be preferentially distributed among a large number 
of HII regions or centrally concentrated toward ring-like 
structures of multiple knots and star clusters. The fact 
that this atypical configuration is found in a galaxy ini- 
tially selected as the host of a rare event such as a GRB 
raises therefore the question of a possible link between 
the physical properties of the WR region and the trigger 
of the long Gamma-Ray Bursts in star-forming galaxies. 

Indeed, the properties observed in this environment 
look similar to a number of specific characteristics com- 
monly found in GRB host galaxies at cosmological dis- 
tances. First, the mid-IR spectroscopic features and the 
SSFR of the WR region both confirm the presence of 
extremely young star formation (i.e., <5Myr) in this 
area. Assuming a consta nt st ar formation rate, the SSFR 
that we derived in Sect. 4.2 indicates that this complex 



of massive star formation would double its stellar mass 
in only 47Myr. This is extremely short compared to 
the typical time scale of galaxy evolution, and although 
it is not as extreme in other GRB host galaxies it is 
consistent with the higher SSFR displayed by this pop- 
ulation when co mpared to field star-forming galaxies at 



similar redshifts (Christensen et al. 20041 Thone et al. 



20081). Second, the WR region has sub-solar metallicity, 
only modest e xtinction by dust and therefore very blue 
optical colors (|Fynbo et aL 2000; Ham mer et al. | |2006| 



Christensen et al.||2 008; Michaiowski et al. 2009J similar 
to what i s typically observed tor GRB hosts at higher 



redsh ifts ( |Fruchter|i9"99}|Le Floc'h et al.|2003{|Han et al 
2010[ ). Among the numerous Hll environments detected 
in the GRB 980425 host it also harbors the largest surface 
brightness region of the galaxy, w hich is consistent with 
the link that was established by Fruchter et al. (2006) 
between the location of GRBs and the very brightest 
regions of the hosts. Finally, optical emission line ra- 
tios like [OIII]/H/3, [Nil] /Ha, [SII]/Ha and [Nil] /[Oil] 
suggest that the ionizing conditions and the stellar pop- 
ulations in the WR region are very similar to that ob- 
served in more distant an d unresolved GRB host galaxies 
( Christensen et aLj2 008 ) . Without claiming a direct link 
between the birth ot the GRB 980425 progenitor and the 
WR region, it would seem therefore difficult from all of 
these observations to argue against a kind of causal link 
between the two. 

Besides, and contrary to the characteristics of the 
WR region, it is worth noting that the integrated prop- 
erties of the GRB 980425 host galaxy differ from those 
of other GRB hosts at high redshift. For example, the 
host of GRB 980425 is an already- evolved galaxy that is 
still forming stars, but at a very low rate wit h respect 
to the stellar mass that it has alrea dy built (Savaglio 
et al ||2009| |Michalowski et al ||2009[ ). Furthermore the 
[OIllj/H/y line ratio measured in the other HII regions of 
the GRB 980425 host are substantially smaller than what 
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is t ypically observed in more distant LGRB host galax- 



ies (Christcnsen et al. 2008L These differences further 



emphasize the connection that must exist between the 
episode of star formation currently taking place in the 
WR region and the trigger of the GRB 980425 in this 
galaxy. 

Nonetheless, the most unsettling fact is that 
GRB 980425 did not occur within the WR region but 
at a projected distance of ~ 900 pc away. What can 
then be the physical connection between the two ? Long 
GRBs represent extremely rare phenomena which are 
more likely to occur from large populations of massive 
stars. Given the very low density of the latter in the 
area where GRB 980425 was found, and considering the 
presence of numerous Wolf-Raye t and OB stars in the 
WR region, Hammer et al. ( |2006[ ) suggested a "run-away 
scenario" in which the GRB 980425 progenitor may have 
come from the WR region itself. This progenitor would 
have been dynamically ejected after a kick or close stellar 
encounters, and traveled at velocities ~ 200-300 km s _1 
across the ISM of the galaxy before ending its life as 
a GRB. Very massive runaway stars with such veloc- 
ities are actually expected from numerical sim ulations 
of dynamical encounters in stellar clusters (e.g., Gvara- 
madze fc Gualandris||2011 ). In the mid-infrared, the de- 
tection of very high-excitation transition lines such as 
[OI V] at 25.9 /jm can indicate the presence o f very hot 
Wolf-Rayet stars ( Schaerer fc Stasihska|[l999| or recent 
supernovae, which in the case of the WR region could 
further strengthen this possible connection with the pro- 
genitor of GRB 980425. In the IRS spectrum however, 
the emission from [OIV] is only observed at the 2a level, 
making this association uncertain. 

Conversely, star formation is still on-going at the lo- 
cation where GRB 98425 was observed and therefore one 
can not fully rule out that its progenitor was truly born in 
this area. For instance, this site could have experienced 
several million years ago a massive starbursting episode 
strongly biased toward the high mass end of the IMF, 
hence explaining the relatively low density of stars now 
observed in this area compared to that of the WR re- 
gion. In this second scenario, it is the physical mech- 
anism that triggered star formation at both locations, 
as well as similarities between the chemical properties of 
these two environments, that would explain a possible re- 
lationship between the WR region and GRB 980425. For 
instance, a single density wave propagating through the 
disc of the galaxy could have triggered over a short time 
scale but in a sequential way the formation of new stellar 
populations along the south-eastern spiral arm, and we 
would be seeing the host while the star-forming episode 
at the GRB location has just terminated. In fact internal 
processes must have played a critical role in driving the 
recent star-forming activity in this object. The galaxy 
ESQ 184-G82 appe ars as an isolated system in the sky 
(Foley et al. 2006) and a trigger of star formation due 
to the accretion of external material via minor merging 
seems quite unlikely. 

6.3. Perspectives 

The current data are unfortunately not sufficient to de- 
termine the exact origin of the progenitor of GRB 980425. 
Further progress may be expected from IFU kinematic 
studies, which could provide better constraints on the 



velocity fields and the dynamical properties of the stellar 
populations in the host galaxy. High-resolution maps of 
the molecular gas distribution inferred from facilities like 
ALMA would also help constraining how the star forma- 
tion efficiency varies in the surrounding of the WR re- 
gion. 

Similarly, characterizing with more statistically rep- 
resentative samples the locations of GRBs with respect 
to the distribution of star-forming regions in their hosts 
should provide more insights into the connection between 
the occurrence of LGRBs and the dynamical ejection of 
massive stars in s tellar clusters. A clear trend has already 
been observed by Fruchter et al. ( 2006 ) between the loca- 
tions of cosmological GRBs ana the brigh test regions of 



their hosts (see also Svensson et al.| [2010). However the 



spatial resolution of current facilities has not yet been 
sufficient to constrain possible offsets at sub-kpc scales. 
Further improvements should thus be expected from next 
generation instruments like JW ST and the Extremely 
Large Telescopes. 

Independently, the properties of the WR region in the 
host of GRB 980425 reinforce the idea that LGRBs are 
triggered within or in the vicinity of extremely young but 
active regions of massive star for mation, lending further 
suppo rt to their collapsar origin (MacFadyen & Woosley 
19991. Similar to what can be int erred from the chem 
ical abundance of the W R region (iHammer et al.||2006 



Christensen et aL]|2008[) this star-formmg activity tends 



to be associated with metal-poor environments, hence 
explaining that LGRBs rather occur within relatively un- 
obscured blue gal axies with sub-solar oxygen abundances 
( |Han et al.|[2010 |. 



7. SUMMARY AND CONCLUSIONS 

Thanks to its rel atively small distance (i.e., 36Mpc, 
Galama et aT||1998| ) the host of GRB 980425 represents 
one of the very tew cases of LGRB hosts where spec- 
troscopic and bolometric properties can be spatially re- 
solved. It thus offers a unique opportunity to con- 
strain the characteristics of the close environment where 
LGRBs take place. In this work we present Spitzer 5— 
160 ^m spectro-imaging data of the GRB 980425 host 
galaxy and the "Wolf-Rayet region" close to which the 
GRB was observed. Our main results can be summarized 
as follows: 

• Given the compactness of the WR region (< 100 pc 
in diameter), the highly- ionized emission lines in 
our IRS spectrum (e.g., [Nelll], [SIV], [Arlll]) sug- 
gest that this environment is associated with ex- 
tremely young (< 5Myrs) massive star formation. 

• While the whole GRB 980425 host galaxy is pro- 
ducing stars at a very low rate compared to its 
already-formed stellar populations, the WR region 
will double its stellar mass in only 47Myrs assum- 
ing constant SFR. It has a total IR luminosity of 
~ 4.6 x 10 8 L and it represents one of the most 
luminous HII regions identified in the nearby Uni- 
verse. 

• In the WR region the relative strength of the PAH 
features over the continuum is similar to what has 
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been observed in other systems with comparable 
oxygen abundance (Z<0.5Z Q ). The extinction 
measured from the optical depth of the 9.7 /im sili- 
cate feature implies modest obscuration (rg 7 Mm ~ 
0.015). 

• Although the WR region produces less than 5% 
of the B-band emission of the GRB host, it con- 
tributes up to 45±10 % of the total IR luminosity of 
the galaxy. Its IR SED also peaks at much shorter 
wavelengths than the global SED of the host, sug- 
gesting that the characteristic temperature in the 
WR region is about twice that of the average dust 
temperature in the interstellar medium of the sys- 
tem. 

• The ratio between the mid-IR and the £>-band lu- 
minosities in the WR region is characteristic of re- 
cent starbursts and much higher than what is ob- 
served in star-forming galaxies with similar Ltr. 
However it is not as extreme as in very young and 
dust-embedded stellar clusters that harbor much 
stronger extinction than derived from the mid-IR 
spectrum of the GRB host. Given the blue col- 
ors of the WR region at optical wavelengths, this 
suggests that this complex of star formation har- 
bors one or several star clusters that have already 
evolved and partly escaped the original molecular 
cloud where they were born. 

The contribution of the Wolf-Rayet region to the bolo- 
metric output of the GRB 980425 host makes this galaxy 
a rather unique object of the nearby Universe. Consid- 



ering that GRBs represent extremely rare events, our 
analysis supports therefore the existence of a causal link 
between the origin of GRB 980425 and the activity of 
star formation in the WR region. More generally, and 
taking also account of the properties characterizing the 
whole population of LGRB hosts, our results corrobo- 
rate the common idea that the triggering of long GRBs 
is mostly associated with very active star formation in 
chemically-young environments. 
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